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Abstract. In this paper, the null controllability in any positive time T of the first-order equation 
(I) x{t) = e ll} Ax(t) + Bu(i) (\9\ < it/2 fixed) is deduced from the null controllability in some positive 
time L of the second-order equation (2) z(t) = Az(t) + Bv(t). The differential equations (1) and (2) 
are set in a Banach space, B is an admissible unbounded control operator, and A is a generator of 
cosine operator function. 

The control transmutation method explicits the input function u of (1) in terms of the input 
function v of (2): u(t,x) = J R k(t, s)v(s) ds, where the compactly supported kernel k depends on 
T and L only. It proves that the norm of a u steering the system (1) from an initial state xo to 
zero grows at most like \\xq\\ exp(a*L 2 /T) as the control time T tends to zero. (The rate a t is 
characterized independently by a one-dimensional controllability problem.) 

In the applications to the cost of fast controls for the heat equation, L is the length of the longest 
ray of geometric optics which does not intersect the control region. 

Key words. Controllability, fast controls, control cost, transmutation, cosine operator function, 
heat equation. 
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1. Introduction. This paper concerns the relationship between the null-con- 
trollability of the following first and second order controllable systems: 



where x and z are the systems trajectories in the Banach space X, xq and zq axe 
initial states, u and v are input functions with values in the Banach space U, A is 
an unbounded generator, B is an unbounded control operator, 9 is a given angle in 
] — 7r/2, 7t/2[, each dot denotes a derivative with respect to the time t and R+ = [0, oo) 
(the detailed setting is given in SJ2J). 

Equation with u = describes an irreversible system (always smoothing) 

and we think of it as a parabolic distributed system with infinite propagation speed. 
Equation l|1.2[l with v — describes a reversible system (e.g. conservative) and we 
think of it as a hyperbolic distributed system with finite propagation speed. For 
example, if A is the negative Laplacian on a Euclidean region and the input function 
is a locally distributed boundary value set by B, then (|1.2jl is a boundary controlled 
scalar wave equation and with 8 = is a boundary controlled heat equation (Sj^l 
elaborates on this example). 

This paper presents the control transmutation method (cf. |5] for a survey on 
transmutations in other contexts) which can be seen as a shortcut to RuselPs famous 
harmonic analysis method in |14j . It consists in explicitly constructing controls u 
in any time T for the heat-like equation in terms of controls v in time L for 

the corresponding wave-like equation H1.2(l . i.e. u(t) = f R k(t, s)v(s) ds, where the 
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x(t) = e w Ax(t) + Bu(t) 
z(t) = Az(t) +Bv(t) 



(teK+), x(Q)=x , 
(t G M), z(Q)=zo, «(0) 



(1.1) 
(1.2) 
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compactly supported kernel k depends on T and L. It proves that the exact control- 
lability of (|1.2|l in some time L implies the null controllability of Ijl.ljl in any time 
with a relevant upper bound on the cost of fast controls for (|l.ip . Thanks to the 
geodesic condition of Bardos-Lebeau-Rauch (cf. for the controllability of the wave 
equation, the application of this method to the boundary controllability of the heat 
equation (cf. JJHJ yields new geometric bounds on the cost of fast controls (extending 
the results of ^I] on internal controllability) . The companion paper J2j concerns the 
quite different case \9\ = ir/2 (in particular there is no smoothing effect). 

The relationship between first and second order controllable systems has been 
investigated in previous papers, always with 9 = and the additional initial data 
i(0) = z\ in X (i.e. considering trajectories of (|1.2(l in the state space X x X). 
In 0, it is proved that the approximate controllability of (11.2(1 with i(0) = Z\ in 
some time implies the approximate controllability of 1(1. lj) for any time (the control 
transmutation method yields an alternative proof), and proves the converse under 
some assumptions on the spectrum of A. The converse is investigated further in |20j 
(in Hilbert spaces) and ^Hj- In a restricted setting, the null controllability of (|l.l[l 
was deduced from the exact controllability of l|1.2Jl with i(0) = z\ in ^1] and ^H] by 
the indirect method of bi-orthogonal bases. 

The study of the cost of fast controls was initiated by Seidman in 16 with a re- 
sult on the heat equation obtained by Russell's method (cf. [S], for improvements 
and other references). Seidman also obtained results on the Schrodinger equation 
by working directly on the corresponding window problem for series of complex ex- 
ponentials (see 52 f° r improvements and references). With collaborators, he later 
treated the case of finite dimensional linear systems (cf. ^S]) and generalized the 
window problem to a larger class of complex exponentials (cf. |17p. The control 
transmutation method generalizes upper bounds on the cost of fast controls from the 
one-dimensional setting (which reduces to a window problem) to the general setting 
which we specify in the next section. 

2. The setting. We assume that A is the generator of a strongly continuous 
cosine operator function Cos (i.e. the second-order Cauchy problem for z(t) = A{t)z is 
well posed and Cos is its propagator). For a textbook presentation of cosine operator 
functions, we refer to chap. 2 of [H] or §3.14 of pQ. The associated sine operator 
function is Sin(i) = J Q Cos(s)ds (with the usual Bochner integral). Cos and Sin are 
strongly continuous functions on R of bounded operators on X. Moreover A generates 
a holomorphic semigroup T of angle ir/2 (cf. th. 3.14.17 of 1 ). In particular S(t) = 
T{e t6 t) defines a strongly continuous semigroup (S(i))tgR + of bounded operators on 
X. In this setting, for any source term / G Lj^^R, X), for any initial data xq, zq and 
z\ in X , the inhomogeneous first and second order Cauchy problems 

x{t) = e l6 A(t)x + f(t) (*GK+), x(0) = z , (2.1) 

z(t) = A(t)z + f{t) (t G R), z(0) = z , z(0) = zi, (2.2) 

have unique mild solutions x G C°(R+, X) and z G C°(R, X) defined by: 

x{t) = S(t)x + [ S(t- s)f(s)dt, z(t) = Cos(t)z Q + Sm(t)z 1 + [ Sm(t - s)f(s)dt. 
Jo Jo 

Remark 2.1. When A is a negative self-adjoint unbounded operator on a Hilbert 
space, T , Cos and Sin are simply defined by the functional calculus as T(t) = exp(iA), 
Cos(i) = cos(ty/~A) and Sin(t) = (V^/T 1 sin^v^)- 
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Following |2I], we now make natural assumptions on B for any initial data in the 
state space X to define a unique continuous trajectory of each system 11.1(1 and (|f .2JI . 
Let A"_i be the completion of X with respect to the norm ||a;||-i = \\(A — [3)~ 1 x\\ 
for some (3 G C outside the the spectrum of A. X-i is also the dual of the space 
Xi defined as the domain D(A) with the norm ||x||i = || (A — (3)x\\. We assume that 
B G C(U, X-i) is an admissible unbounded control operator hi the following sense: 



Vt > 0,Vu G L 2 ([0,t}-U), / S(s)Bu(s)ds G X and / Sin(s)Bit(s)ds G X, (2.3) 



In this setting, for any xq and zq in X, for any u and i> in L 2 oc (R + ; U), the unique 
solutions a; and z in C°(]R; A") of (f 1 . 1|> and (|1.2|l respectively are defined by: 



The natural notions of controllability cost for the linear systems and (|1.2[) are: 

Definition 2.2. The system $1.1)) is null-controllable in time T if for all xo in 
X, there is a u in L 2 (M.;U) such that u{t) = for t ^ [0,T] and x{T) = 0. The 
controllability cost for $1.1)) in time T is the smallest positive constant k± : t in the 
following inequality for all such <j)Q and u: ^\\u{t)\\ 2 dt ^ ki,t||^o|| 2 - 

The system M.2\) is null-controllable in time T if for all zq in X , there is a v in 
L 2 (K; U) such that v(t) = for t <£ [0,T] and z(T) = z(T) = 0. The controllability 
cost for $!.!&) in time T is the smallest positive constant K2.T in the following inequality 
for all such zq and v: J^\\v(t)\\ 2 dt ^ K2,t||2o|| 2 . 

Remark 2.3. Equivalently, for all xt in S(T)X, there is a u in L 2 (0,T; U) such 
that x(Q) = and x(T) = xt, and, for all zq and zt in X , there is a v in £ 2 (0, T; U) 
such that i(0) = z{T) = 0, z(0) = zq and z(T) = zt- 

3. The results and the method. Our estimate of the cost of fast controls for 
Hl.l|) builds, through the control transmutation method, on the same estimate for 
a simple system of type <f 1 . 1 1) . i.e. on a segment [0,L] with Dirichlet (N — 0) or 
Neumann (N = 1) condition at the left end controlled at the right end through a 
Dirichlet condition: 

M = e l6 d 2 <P on ]0,T[ x ]0,L[, d? ^ s=0 = 0, ^ s=z , = u, lt=o - fa. (3.1) 

With the notations of x = tj>, A = d 2 s on X = L 2 (0,L) with D(A) = {/ G 
H 2 (0, L) | <9f /(0) = f(L) = 0}, || - 1| i with (3 = is the homogeneous Sobolev ij 2 (0, L) 
norm, and B on U = C is the dual of C G C(X%; U) defined by Cf = d s f(L). 

It is well-known that the controllability of this system reduces by spectral analysis 
to classical results on nonharmonic Fourier series. The following upper bound for the 
cost of fast controls, proved in is an application of a refined result of Avdonin- 
Ivanov-Seidman in |17| . 

Theorem 3.1. There are positive constants a and 7 such that, for all N G {0, 1} ; 
L > 0, T G ] 0, inf (7r, i) 2 ] , the controllability cost kl.t of the system $3.1]) satisfies: 
kl,t ^ 7exp(ai 2 /T). 

This theorem leads to a definition of the optimal fast control cost rate for (13.1(1 : 

Definition 3.2. The rate a* is the smallest positive constant such that for all 
a > a» there exists 7 > satisfying the property stated in theorem V3.1\ 
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Remark 3.3. Computing a* is an interesting open problem and its solution does 
not have to rely on the analysis of series of complex exponentials. The best estimate 
so far is a* G [1/2, 4(36/37) 2 ] for 6 = (cf. \n$). 

Our main result is a generalization of theorem 13.11 to the first-order system (ll.lfl 
under some condition on the second-order system l|1.2[) : 

Theorem 3.4. If the system \l.ty is null- controllable for times greater than L*, 
then the system 11. 1\) is null- controllable in any time T . Moreover, the controllability 
cost Ki t T of l\l.l\) satisfies the following upper bound (with a* defined above): 

limsupThiKi^ ^ (3-2) 

Remark 3.5. The upper bound \°3.Hfy means that the norm of an input func- 
tion u steering the system l^l.l)) from an initial state xq to zero grows at most like 
7||xo|| exp(aL 2 / (2T)) as the control time T tends to zero (for any a* and some 7 > OJ. 
The falsity of the converse of the first statement in th. \3.4\ is well-known, e.g. in the 
more specific setting of 

Remark 3.6. As observed in J^, \3.^ yields a logarithmic modulus of continuity 
for the minimal time function T m i n : X — > [0, +00) of l\l.l\) : i.e. T m i n (xo), defined 
as the infimum of the times T > for which there is a u in L 2 (R;Y) such that 

Jo II^WII ^ 1> U W = f° r t 4- [0;^1 an d X {T) = 0, satisfies: for all a > a*, 
there is a c > such that, for all xg and x' in X with \\xq — x' \\ small enough, 
\T min (x ) - T min (x' )\ < all/ ln(c/||x - x' \\). 

It is well known that the semigroup T can be expressed as an integral over the 
cosine operator function Cos (cf. the second proof of th. 3.14.17 in 0): 

VxEXyteC s.t. |arg*| < tt/2, T{t)x = J k(t, s) Cos{s)xds , (3.3) 

where fc(0, s) = 5(s) and k(t, s) = exp(— s 2 / (At))/ \/irt for Ret > 0. This equation has 
been referred to as the abstract Poisson or Weierstrass formula. Starting with the 
observation that k is the fundamental solution of the heat equation on the line, i.e. k 
is the solution of dtk = d^k with the Dirac measure at the origin as initial condition, 
the transmutation control method consists in replacing the kernel k in Ij3.3|l by some 
fundamental controlled solution on the segment [-L, L] controlled at both ends (cf. 
(I5.8[) h The one dimensional th. 13.11 is used to construct this fundamental controlled 
solution in 21 and the transmutation is performed in f|5l 

4. The fundamental controlled solution. This section begins with an outline 
of the standard application of ^2] to the proof of th. 13.11 Following closely §5 of jll| . 
the rest of the section outlines the construction of a "fundamental controlled solution" 
k in the following sense, where T>'(0) denotes the space of distributions on the open set 
O endowed with the weak topology, Ai(0) denotes the subspace of Radon measures 
on O, and S denotes the Dirac measure at the origin: 

Definition 4.1. The distribution k e C°([0,T];M(] — L,L[)) is a fundamental 
controlled solution for 14- 1\) at cost (7,0;) if 

d t k = e- w d 2 s k mVQQ,T[x]-L,L[) , (4.1) 
fc lt=0 =<5 and fc lt=T = 0, (4.2) 

IWIiw[x]-L,L[) ^7e Qi2/T (4.3) 
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The operator A defined at the beginning of is negative self-adjoint on the 
Hilbert space L 2 (0, L). It has a sequence {^n}neN* of negative decreasing eigenvalues 
and an orthonormal Hilbert basis {e n } ne N* in L 2 (0,L) of corresponding eigenfunc- 
tions. Explicitly: y/—JI^ = (n + v) tt/L with v = for N — (Dirichlet) and v = 1/2 
for AT = 1 (Neumann). First note that th. l3.H can be reduced to the case L = n by the 
rescaling (t, s) h- > (a 2 t,as) with a = L/71". In terms of the coordinates c = (c„)„ e pf* 
of A N / 2 fo in the Hilbert basis (e n )„ 6 N* where fo is the initial state of the dual ob- 
servability problem, th. 13.11 with L = 7r reduces by duality to the following window 
problem: 3a > 0, 3 7 > 0, VT G ] 0, tt 2 ] , 



Vc G i 2 (N*), ^ |c„| 2 ^ 7 e Q7r2/T / |F(*)| 2 cft where F(t) = ^ c, 

nSN* ^° n=l 



e exp(i6)/i„ 



Since this results from th. 1 of [TJ| with A„ ~ ie l8 n 2 as in §5:2 of the proof of 
th. 13. II is completed. 

Now we consider a system governed by the same equation as 1(3. lj) but on the 
twofold segment [—L,L] controlled at both ends: 

d t( f>-e- ie d 2 ( f> = Q m]0,T[x]-L,L[, <^ 3=±L = u±, lt=o = <fo, (4.4) 

with initial state (fio G £ 2 (0, L), input functions U- and w+ in L 2 (0, T). As in propo- 
sition 5.1 of |11|. applying th. l3.1l with A^ = to the odd part of <po and with AT = 1 to 
the even part of </>o proves that the controllability cost of (|4.4(l is not greater than the 
controllability cost of 1(3.1(1 and therefore satisfies the same estimate stated in th. 13.11 
As in proposition 5.2 of |11| . we may now combine successively the smoothing effect of 
1(4.4(1 with no input (i.e. u + = u_ = 0) and this controllability cost estimate (plugged 
into the integral formula expressing <fi in terms of 0o and u± = 4>] s= ±l) to obtain: 

Proposition 4.2. For all a > a*, there exists 7 > such that for all L > and 
T G ]0, inf (n/2, L) 2 ] there is a fundamental controlled solution for j4-4\ ) at cost (7, a) 
(cf. def.g^. 

5. The transmutation of second-order controls into first-order controls. 

In this section we prove th. 13.41 

Let xq G X be an initial state for i|l.l|) and let L > L*. Let z G C°(R+; X) and 
v G L 2 (K + ;U) be the solution and input function obtained by applying the exact 
controllability of p. 2(1 in time L to the initial state zq — xq. 

We define z_ G C°(R;X) and uG L 2 (R ; [/) as the extensions of £ and v by reflection 
with respect to s — 0, i.e. z(s) = £( s ) = z(—s) and v(s) — v(s) — v(—s) for s ^ 0. 
They inherit from l|2.4|) : 

z(t) = Cos(i)x + / Sin(t - s)Bv(s)ds. (5.1) 
Jo 

Def. 12.21 of K2,l implies the following cost estimate for v; 

\\v(s)\\ 2 ds = 2 [ \\v(s)\\ 2 ds < 2 K2 , i ||x || 2 . (5.2) 



Since Z) (A) is dense in A", there is a sequence (x n )n£N* in -D(A) converging to ceo in 
X. Since A"i is dense in X-i, there is a sequence (/ n ) n eN* m C 1 (M; Xi) converging 
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to By in L 2 (R; For each neH', let z n be denned in C 2 (R; X) by: 

2 n (t) = Cos(i)x„ + / S'm(t - s)f n (s)ds, 



which converges to z(t) in X for all t due to 1)5. Since z n is a genuine solution of 
z(t) = Az n {t) + /„(*) (cf. lem. 4.1 of [6]), we have for all 93 in D(A'): 

d 2 

(2»00,¥>) G # ( R ) and 4^ ><P) A' (p) + {f n (s),<p). 

Hence, (z n (t),(p) = (x n ,tp)+ / (f-s)(z n (s),iV)+ / ¥>)• Passingtothc 

Jo Jo 

limit, yields (z(t), tp) = (xq, tp) + / (t-s)(z(s),iV)+ / s )(By(s), ip) . Therefore: 

Jo Jo 

s h-> (z( s ), ip) e H 2 (R) and — («(«), = (z(s),Aip) + (By(s),<p), (5.3) 

(z(s),ip)=0 and 4fe( s ):^>=0 for Isl = L. (5.4) 
as 

Let a > a* and T g]0, inf (1, L 2 )[. Let 7 > and fc G C°([Q,T];MQ - L,L[)) 
be the corresponding constant and fundamental controlled solution given by propo- 
sition We define k G C°(R+; A-l(K)) as the extension of k by zero, i.e. k(t, s) = 
k(t,s) on [0, T] x] — L, i[ and k is zero everywhere else. It inherits from k the following 
properties 

d t k = e l9 d 2 k mV'(]0,T[x] -L,L[) , (5.5) 
k]t=o = ^ and &t=T = > (5-6) 
Hfc|lW R) <7 e QL2/T (5-7) 

The main idea of the proof is to use k as a kernel to transmute z and y into a 
solution x and a control u for 1)1.1)) . The transmutation formulas: 

= J k(t,s)z(s)ds and Vt > 0, «(*) = J k(t, s)y(s) ds, (5.8) 

define x G C°(M+;X) and u £ L 2 (R+;U) since k G C°(R + ; 7W(R)) n L 2 (M+; L 2 (M)), 
z G C°(R;X) and w G L 2 (R;t/). The property of fc implies x(0) = x and 

x(T) = 0. Equations 1)5.3)1 . 1)5.4)1 and 1)5. 5|l imply, by integrating by parts: 

Vtp G £>(^4'), * ^ ¥>) G and 4 M*), ¥>) = A V) + (-Bu(t), V?). 

This characterizes x as the unique solution of 1)1.1)1 in the weak sense (cf. 0), which 
implies that x and u satisfy 1)2.4)1 . Since J\\u(t)\\ 2 dt < JJ\k(t, s)\ 2 dsdt J\\y(s)\\ 2 ds, 
1)5.7)1 and 1)5.2)1 imply the following cost estimate which completes the proof of th. 13.41 

T \\u{t)\\ 2 dt<:2^ Ll e aL2 ' T \\x \\ 2 . 
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6. Geometric bounds on the cost of fast boundary controls for the heat 
equation. When the second-order equation ljl.2jl has a finite propagation speed and 
is controllable, the control transmutation method yields geometric upper bounds on 
the cost of fast controls for the first-order equation . From this point of view, this 
method is an adaptation of the kernel estimates method of Cheeger-Gromov- Taylor 
in 0. This was illustrated in JI] and on the internal controllability of heat 
and Schrodinger equations on Riemannian manifolds which have the wave equation 
as corresponding second-order equation. Some similar lower bounds are proved in 
these papers (without assuming the controllability of the wave equation) which imply 
that the upper bounds are optimal with respect to time dependence. In this section, 
we illustrate the control transmutation method on the analogous boundary control 
problems for the heat equation. 

Let (M, g) be a smooth connected compact rt-dimensional Riemannian manifold 
with metric g and smooth boundary dM . When dM ^ 0, M denotes the interior 
and M = M U dM. Let A denote the (negative) Laplacian on (M, g) and d v denote 
the exterior Neumann vector field on dM. The characteristic function of a set S is 
denoted by \s- 

Let X = L 2 (M). Let A be defined by Af = Af on D(A) = H 2 (M) n H%(M). 
Let C be defined from D{A) to U = L 2 (dM) by Cf — <9i,/]r where F is an open 
subset of dM, and let B be the dual of C. With this setting, ljl.lt with 6 = and 
<jl.2jl are the heat and wave equations controlled by the Dirichlet boundary condition 
on r. In particular ljl.2jl writes: 

d 2 z -Az = on I f xl, z = \tv on R t x dM, 
z(0)=z eL 2 (M), C(0)=0, veL 2 oc (R;L 2 (dM)), ( "''' 

It is well known that B is an admissible observation operator (cf. cor. 3.9 in 3 ). To 
ensure existence of a null-control for the wave equation we use the geometrical optics 
condition of Bardos-Lebeau-Rauch (specifically example 1 after cor. 4.10 in [3]): 

There is a positive constant Ly such that every generalized geodesic of length , . 
greater than Ly passes through T at a non-diffractive point. 

Generalized geodesies are the rays of geometrical optics (we refer to |l()j for a presen- 
tation of this condition with a discussion of its significance). We make the additional 
assumption that they can be uniquely continued at the boundary dM. As in 0, to 
ensure this, we may assume either that dM has no contacts of infinite order with 
its tangents (e.g. dM = 0), or that g and dM are real analytic. For instance, we 
recall that (j6.2|) holds when T contains a closed hemisphere of a Euclidean ball M of 
diameter Lp/2, or when T = dM and M is a strictly convex bounded Euclidean set 
which does not contain any segment of length Ly. 

Theorem 6.1 ( 3 ). If holds then the wave equation |6'.i|) is null- controllable 
in any time greater than Ly- 

Thanks to this theorem, th. 13.41 implies: 

Theorem 6.2. If holds then the equation: 

d t x - e ie Ax = on R t x M, x = xm on R t x dM, 
x(0) =x Q £ H- 1 {M), u G Lf oc (R;L 2 (dM)), 

is null- controllable in any time T . Moreover, the controllability cost k\_t (cf. def. \2.ty) 
satisfies (with a* as in def. VS.ty) : limsupThiK^T ^ a*Lp. 
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